Resonant Fibonacci Quantum Well Structures 
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We propose a resonant one-dimensional quasicrystal, namely, a multiple quantum well (MQW) 
structure satisfying the Fibonacci-chain rule with the golden ratio between the long and short inter- 
well distances. The resonant Bragg condition is generalized from the periodic to Fibonacci MQWs. 
A dispersion equation for exciton-polaritons is derived in the two-wave approximation, the effective 
allowed and forbidden bands are found. The reflection spectra from the proposed structures are 
calculated as a function of the well number and detuning from the Bragg condition. 

PACS numbers: 42.70.Qs, 61.44.Br, 71.35.-y 



INTRODUCTION 

The concept of quasicrystal as a non-periodic struc- 
ture with perfect long-ranged bond orientational order 
was brought in solid-state physics by Levine and Stein- 
hardt [1|. It was cxtcndcd to optics in Ref. |2|, whcrc 
a one-dimensional (ID) quasicrystal model constructed 
of dielectric layers forming the Fibonacci sequence was 
proposed. At just the same time the concept of photonic 
crystals was suggested by Yablonovich [3| and John [j]. 
Since then the ID photonic Fibonacci quasicrystals have 
been extensively studied[5|, la, |7|, |8[ . 

In this paper we introduce a new nanoobject, the Fi- 
bonacci quantum well (QW) structure with inter-well 
spacings arranged in the Fibonacci sequence. This means 
that the thickncss of barriers separating the wcUs can 
take one of two values so that the ratio between the 
long and short inter-well spacings cquals the golden mean 
T = (V5+l)/2. We focus on the Hght propagation in such 
a medium in the frequency region around the resonance 
frequency ujq of a two-dimensional exciton in the quan- 
tum well. The barriers are assumed to be thick enough so 
that the excitons in diffcrcnt wells are coupled only via 
electromagnetic field. Thus, the object under study is 
a resonant photonic quasicrystal, an intermediate struc- 
ture between completely ordered and disordcred media, 
namely, periodic MQWs with a fixed inter-well spacing 
and MQWs with random inter-well spacing. 

Among periodic QW structures, of particular interest 
are the resonant Bragg structures with the period satis- 
fying the Bragg condition 



q{uja)d = TTJ, j = 1,2. 



(1) 



where q{uj) = Lorih/c and q{too) is the light wave vector 
at the exciton resonance frequency loq, Uf, is the back- 
ground refractive index of both QW and barrier mate- 
rials, d is the structure period, and e is the light ve- 
locity. The periodic resonant BraggMQWs have been 
first considered theoretically in Ref. [9| and then investi- 
gated in a number of theoretical as well as experimental 
Works [iS im, [il [il Ei, m, [il m . It was established 



that, for small enough numbers N of QWs (superradiant 
regime), the optical reflection spectrum is described by 
a Lorentzian with the halfwidth NTq + T, where Tq and 
r are, respcctively, the exciton radiative and nonradia- 
tive damping rates in a single QW [9|, [ij]. For a large 
number of wells (photonic crystal regime) , the reflection 
coefficient is dose to unity within the forbidden gap for 
exciton polaritons propagating in infinite periodic system 
and rapidly decreases near the gap edges wn — A A/j and 
A/v7, where A = VWV^ Ìli, Q [Ìl, lÌÌ • 



Wo 



In Section II we will show that a generalized resonant 
Bragg condition analogous to Eq. (JT]) can be formulated 
for the resonant Fibonacci MQW structures, although 
the latter are aperiodic. In Section III the significance of 
the proposed condition is verified by numerical calcula- 
tions of the reflection spectra from the structures tuned 
on and slightly detuned from this condition, and the de- 
pendence of the reflection spectra on the number of wells 
is analyzed and compared with those for the periodic 
Bragg structures. In Section IV we apply a two-wave ap- 
proximation in order to determine the band gaps in the 
exciton-polariton spectrum of the Fibonacci structures 
and show that the simple analytic theory allows one to 
interpret quite well the numerical results. 
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FIG. 1: Scheme of the Fibonacci QW structure jFg with A'' = 9 
QWs. 



RESONANT BRAGG CONDITION FOR 
FIBONACCI MQW STRUCTURE 

The structure under consideration is schematically de- 
picted in Fig. 1. It consists of iV identical QWs embedded 
in a matrix with dielectric Constant £{,. The inter-weU dis- 
tances take two values represented by long and short seg- 
ments of length l and s, respectively. For the Fibonacci 
chain, the coordinate Zm of the m-th QW center is givcn 
byH 



Zm = d{m- 1) + 



7-{t})- 



(2) 



where the integer m runs froni 1 to iV, r is the golden 
ratio, d is the averagc pcriod of the structure givcn by the 
product s(3 — r), and {x} is the fractional part of x. An 
alternative way of deflning Zm is bascd on the rccurrencc 
relation J-j+i = {J^j,!Fj-i} for finite Fibonacci chains of 
the order j + 1 , j and j — \, with initial conditions J-i ~ S, 
Ti = L, where S and L are the segments with lengths s 
and l = Ts, respectively. [18[ Then z„i are coordinates of 
boundaries between the segments in the Tj sequence. 

The exact reflection coefficicnt of the light normally 
incident on such a structure froni the left half-space can 
be obtained by standard transfer matrix method [l9| . In 
order to forni the base for formulation of the resonant 
Bragg condition for the Fibonacci structurcs wc will an- 
alyze the reflection in the first-order Born approximation 
neglecting multireflection processes and summing up the 
amplitudes of waves reflected from distinct wcUs. Then 
the amplitude reflection coefficient, rpf{uj), from the A^- 
well Fibonacci structure at the light frequency co is given 

by 

rjvM«iV/[gH,Af]riH, (3) 

whcrc /(q, N) is the structure factor of the system, 

N 

N 



1 ^ 



and ri is the reflection coefficient from a single QW, 

irò 



ri{u)) = 



u!o- uj -i(r + ro) 



For the semiinfinite Fibonacci MQWs the structure fac- 
tor f{q) = limjv^oo fiOjN) can be presented in the fol- 
lowing analytical forni [l8| 



/(g) 



7 , S2q,G^,^,fhh', Ghh' 
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{h + h'/r), 



h,h' — — c 



sin Shh' /.r-2 
Jhh' = —5 exp 1 Shh' 



Shh' 



(4) 
(5) 



S, 



hh' 



r^ + 1 



{tìi' -h). 



Allowcd diffraction vectors Ghh' form a dense pseudo- 
continuous set. The largest values of \fhh'\ are reached 
for the pairs {h,h') coinciding with two successive Fi- 
bonacci numbers (-F,, Fj^i) with Fj dcfincd rccursively 

Thus, for 



F. 



j-i- 



by i^o = 0, Fi = 1 and Fj+i = Fj 
[Kh') = {F,,F,^,) = (1,0), (1,1), (2,1), (3, 2) and (5,3) 
corresponding to j ~ 1...5, the modulus of fhh' equals to 
« 0.70, 0.88, 0.95, 0.98 and 0.99, respectively. For {h, h') 
iiot belonging to this particular set, values of \,fhh'\ are 
signiflcantly siiialler. It foUows then that if the exciton 
rcsonancc frequency satisfics the condition 



Lo^nb 



d = Tri F, 



F,-i 



.7 = 1,2.. 



(6) 



the coefficient ([3]) at w = o^q and large N amounts to 



rN = Nf[q{uo),N]ri{uo) 



NTof, 



hh' 



Fo + r 



This is of the same order of magnitudo as the reflection 
coefficient calculated in the same Born approximation for 
a pcriodic resonant Bragg structure satisfying Eq. ([T|). 
Hence Eq. ^ is indccd a resonant Bragg condition gcn- 
cralized for the Fibonacci MQWs. In the foUowing we fix 
the value of loq, consider the average period d as a vari- 
ablc paramctcr and use the notation dj for d given by 
Eq. ([6]) for the integer j. The corresponding thicknesses 
Sj , Ij of the short and long segments are related with dj 

by 



dj/i3 - t) , Ij = rfjT/(3 - t) 



(7) 



The estimation ([3]) for rjv is valid until |rjv| <C 1, i.e., 
if A^Fo <C maxjlwo — w|,r}. Otherwise one has to take 
info account the multireflection of the light waves from 
QWs which is readily achieved by the standard transfer- 
matrix numerical calculation. The results are presented 
and analyzcd in the ncxt scction. 



CALCULATED REFLECTION SPECTRA 

Figure 2 prescnts reflection spectra calculated for four 
structures containing A'^ = 50 quantum wcUs. The ex- 
citon parameters used are as follows: fuo^ ~ 1.533 cV, 
^Fo = 50 /xeV, hT = 100 /.tcV, m == 3.55. Curve 1 is 
calculated for the resonant Fibonacci QW structure sat- 
isfying the exact Bragg condition ^ with j = 2 so that 
d = di, s ~ S2 and l ~ li- Curves 2 and 3 correspond 
to the Fibonacci structures with the barrier thicknesses 
slightly dctuned from S2 and li'- s/s2 = l/h = 1-02 for 
curve 2 and s/s2 = l/h = 0.98 for curve 3. Curve 4 
describes the reflection from the periodic Bragg struc- 
ture with the same exciton parameters and the period 
d = TT/q{LUo)j satisfying Eq. ([T]). From comparison of 
curves 1 and 4 we conclude that the reflection spec- 
tra from the resonant periodic and Fibonacci structures 
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FIG. 2: Reflection spectra calculated for three Fibonacci 
structures satisfying the resonant Bragg condition (|6]) for 
j — 2 (curve 1) and detuned by ±2% from this condition 
(curves 2 and 3) in comparison with the reflection spectrum 
from the periodic resonant Bragg QW structure (curve 4). 
The values of parameters are indicated in text. 



tuncd to the Bragg conditions ([T]) and © are closc to 
each other outside the frequency region around wq- More- 
over it foUows from curves 2 and 3 that a shght deviation 
from the condition (O results in a radicai decrease of the 
effective spectral halfvifidth. This sensitivity is a charac- 
teristic feature of Fibonacci structures, as it happens for 
the periodic Bragg QW systems. The remarkable struc- 
tured dip in the middle of the spectrum 1 is the only 
qualitative difference from the periodic structures, the 
origin of this dip is explained in the next section. Now 
we turn to analysis of reflection spectra as a function of 
the QW number TV and index j in Eq. ©. 

Evolution of the reflection spectra with the QW num- 
ber N is illustrated in Fig. 3a. The spectral envclopc 
smoothed to ignore dip in the middle shows a behavior 
similar to that of the conventional Bragg QW structure. 
Indccd. for small N the envelope is a Lorentzian with the 
halfwidth incrcasing as a linear function of N. This is 
a straightforward manifestation of superradiant regime, 
which, as we can see here, does not necessarily require pc- 
riodicity even if the inter-well distances are comparable 
to the light wavelength. The saturation of the spectral 
halfwidth (photonic crystal regime) begins at large N of 
the order of y^ujo/TQ, in a similar way as for the periodic 
Bragg structures. The shape of the spectra for large N 
allows US to suppose existence of two wide symmetrical 
stop bands in the energy spectrum of the structure with 
an allowed band between them. Of course, the applica- 
tion of terms "allowed" and "stop" bands to an apcriodic 
structure is questionable. In section IV we show that 
nevertheless these terms are applicable in a reasonable 
approximation. 

Figure 3b prcsents the reflection spectra of the Fi- 



FIG. 3: Reflection spectra from the resonant Fibonacci 
structures. (a) Six curves are calculated for the struc- 
tures satisfying the condition ((6)l with j = 2 and A'^ = 
20, 50, 80, 110, 150, 200. The number of wells is indicated near 
each corresponding curve, (b) Curves 1, 2, 3 are calculated 
for the structures with A^ — 200 and indices j — 2, 3, 5 in 
Eq. ([SJ. Vertical lines connected by a horizontal bar indicate 
the exciton-polariton high-frequency gap given by Eq. (|10|l . 
Other parameters are the same as in Fig. 2. 



bonacci QW structures containing a large number of 
wells, A^ = 200, and satisfying Eq. ^ with three dif- 
ferent values of j. AH the curves indicate an existence of 
the stop and allowed bands. Howevcr, the band widths 
are j-dependent: the stop band (or gap) indicated by a 
unitcd pair of vertical lines and the middle dip are both 
squeezed with increase of j . 



EXCITON-POLARITON ENERGY SPECTRUM 

For the light propagating in a system of identical QWs 
located at the points z„i (m = 1, 2...), the equation for 
the clectric field can be written as [l9Ì 



il 



E{z) 



2qT„ 



LOn 



ir 



^(5(z- Z„i)E{Zm) 



(8) 



whcrc q = q{uj) and we assume that quantum wells are 
thin as compared to the light wavelength. We consider 
a semiinfinite Fibonacci QW structure with the average 
pcriod dj satisfying Eq. ([6|) for a certain vaine of j. Using 
the above mentioned properties of the coefhcients f^w in 
the structure factor f{q) we can retain in the sum (|4]) only 
one term fhh'S2q,G^^, with {h,h') = {Fj,Fj^i). In other 
words we take info account only one diffraction vector 
Ghh' corresponding to the condition ([6]) and neglect ali 
other possiblc diffraction vectors. In this approximation 
we can prcscnt the light wave as a sum of two piane wavcs 



with the wave vectors K and K' ~ K — Ghh' assuming 
K « Ghh' l'^'- The amplitudes of the chosen spatial har- 
monics, Ek and Ek'^ satisfy the following two coupled 
equations 



{q^-K^+X)EK + xfhh'EK'^Q, 
xfhh'Eii + {q^ - K'^ + x)Ek' = , 



(9) 



where 



X = 



2<?ro 



d{ujQ — uj — ir) 



In the following analysis we ignore the exciton dissi- 
pation, ncglccting the nonradiative damping. Thus the 
frcquency axis is divided into intervals of purcly allowed 
and forbidden bands with propagating and evanescent 
polaritonic solutions. In the allowed bands the solutions 
are characterized by real values of the wave vector K. It 
is convenient to reduce the exciton-polariton dispersion 
oj{K) to the "first Brillouin zone" defined in the interval 
-Ghh' 1^2. < K < Ghh'/2. The detailed behavior oi uj{K) 
inside this interval lies out of the scope of the present 
paper. Note that, in dose vicinity to wq, the two-wavc 
approximation is inadequate and the polariton dispersion 
should be calculated taking into account an admixture of 
a lot of piane waves. Here we consider only the exciton- 
polariton eigenfrequencies at the edge of the Brillouin 
zone, K = -K' = Ghh'/2. It follows from Eq. ® that 
four eigenfrequencies at this point are given by 



■ UJQ àz Ai 



I l + \.fhh'\ 

2 {h + h'/r) ' 



(10) 






2(/i + /i'/r) 

In accordance with Fig. 3 we attribute the interval to^ < 
UJ < w^j to the exciton-polariton upper stop band (la- 
belled by index "-I-") and the interval between a;"^^ and 
Wj" to the lower stop band (labelled by index "-"). The 
subscripts "in" and "out" denote the stop-band edges 
inner and outcr with rcspcct to ujq. 

The values of uj^ and w^j are marked by vcrtical lines 
in Fig. 3b. One can see an excellent agreement between 
the band edges revealed in the calculated spectra and 
thosc given by Eq. (|10p which unambiguously confirms 
the interpretation of the frcquencies ([T0|) . 

Equations pO|) can bc rcduced to thosc for the peri- 
odic resonant Bragg structures as soon as \fhh'\ is set 
to unity and h + h' /t is replaced by the integer j. For 
\fhh'\ = 1 the inner eigenfrequencies merge at ujq and a 
single band gap of width 2A/y(J is formed. In the Fi- 
bonacci QW structures \fhh'\ < 1 and, as a result, an 
allowed band opens between cJ;^ and uj^. We note that a 
qualitatively similar band structure can be realized when 
the periodic MQWs has a compound elcmentary celi. 15 1 



One can easily show that, also in this case, the modulus 
of the structure factor is smaller than unity. Morcover, 
Eqs. (UHI) can be reduced to Eqs. (26) of Ref. [H if 
h + h'/r and \fhh'\ are replaced, respectively, by 1/2 and 
I cos qd2 1 , where d2 is the Inter- well distance in the com- 
pound unit celi of a periodic structure with two QWs in 
the supercell. 

In the Fibonacci QW structure the decrease of stop- 
band widths with the increasing h + h'/r is related to the 
corresponding increase of the average period d in Eq. ^ 
and it is analogous to the j~^^^ power law of the band 
width for the periodic resonant Bragg structures. The 
middle allowed band width decreases even faster because, 
as mentioned above, the value of \fhh'\ tends to unity 
and, therefore, the value of -^/l — \fhh'\ rapidly vanishes 
as the index j in Eq. ^ changes from 2 to 5. 



CONCLUSIONS 

We have introduccd into considcration resonant ID 
photonic quasicrystals based on Fibonacci QW struc- 
tures. The analysis of light reflection in the Born approx- 
imation has been used to formulate the resonant Bragg 
condition for this system. The results of straightforward 
transfer-matrix numerical calculation confimi the rele- 
vance of the generalized Bragg condition imposed on the 
aperiodic system under study. For a small number TV 
of QWs, the Fibonacci structures show the superradiant 
behavior while, for high values of N exceeding a/wq/Fo, 
the photonic crystal regime with distinct stop bands in 
optical spectra is reached. A qualitative difference with 
respect to the periodic resonant Bragg QW structures lies 
in the presence of a structured dip in the reflection spec- 
trum around the exciton resonance frequency ujq . An ap- 
proximate two-wave exciton-polariton model allows one 
to describe the widths of the allowed and forbidden bands 
as a function of the structure paramctcrs. 

This work was supported by RBFR and the "Dynasty" 
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